Spin orbit induced torque in collinear spin valve structure and associated entropy 
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We predict that due to spin-orbit(SO) interaction there is a torque even for the parallel configura- 
tion of spin valve structure (F1/2DEG/F2). This torque arises due to spin orbit interaction. We 
develop a scattering theory for spin density matrix which allows us to study this effect quantita- 
tively. Further we show that the von-Neumann entropy associated with transport of polarization 
can show non-linear behavior as a function of absolute angle and oscillator behavior as a function 
of SO interaction strength. 

PACS numbers: 72.25-b,72.25.Dc, 72.25.Mk 
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Since the prediction of spin transfer torque phenomena 
ten year ago by Slonczewski [1] and Berger [2], consider- 
able experimental evidence has been accumulated that a 
spin-polarized current can induce switching or precession 
of the films magnetization [3] [4] [5] [6] [7]. In their pio- 
neering work Slonczewski and Berger predicted for spin 
valve geometry that for uniform and non-collinear mag- 
netization in two ferromagnetic layers , magnetization of 
free layer will precess due to the spin transfer torque. 
This effect vanishes for collinear or parallel magnetiza- 
tion since the spin transfer torque varies as sin(0) where 
9 is the angle between the magnetization of two ferro- 
magnetic layers (see eq. 8 in Ref. [2]). 

On theoretical side a lot of work has been done finding 
new phenomena related to spin transfer torque for e.g. 
dynamic exchange coupling in magnetic bilayers [8] , spin 
echo [9] etc. On the other hand the phenomena has been 
modeled by variety of ways including, combined ballistic 
and diffusive model [10], first principal calculation [12] 
microscopic spin transfer model [13], and magnetization 
dynamics [11]. 

All these studies were done in absence of spin-orbit 
(SO) interaction which is rather important in low dimen- 
sional heterostructure. In presence of SO interaction the 
situation differs dramatically. To elaborate this point we 
specifically consider a two dimensional heterostructure 
consisting of a two dimensional electron gas sandwiched 
between two ferromagnetic contacts (F1/2DEG/F2). In 
such a system with a prominent source of SO interac- 
tion is asymmetric confining potential which gives rise to 
Rashba spin-orbit (RSO) interaction and it magnitude 
can be varied by changing the external electric filed ap- 
plied for confining the structure. We take the plane of 
2DEG as xy plane and a normal to it as z axis which de- 
fines the coordinate system which will be used through- 
out this paper. In this frame the RSO interaction can be 
written as 



H so = B R (k) • er, 



(1) 



where BR,(k)=A(fcj,x — k x y) and <x is a vector of Pauli 
matrices. Magnitude and direction of Br(1c) depends on 



instantaneous momentum vector, infact the direction is 
perpendicular to instantaneous wave vector. 

When an electron is injected into 2DEG its spin pre- 
cesses due to the RSO given by eq.(l) as long as injected 
spin direction is not an eigen state of eq. (1). In hy- 
brid structurer like F1/2DEG/F2 electrons are injected 
not only perpendicular to the interface rather over all 
possible angles. Hence the RSO will cause spin preces- 
sion which implies that when electrons reaches detector 
ferromagnet F2 it spin will develop a component perpen- 
dicular to the magnetization of F2 which is parallel to 
the magnetization of Fl. Hence this perpendicular com- 
ponent will give rise to a torque on the magnetization F2 
and vice versa. This torque arises due to the spin pre- 
cession caused by RSO interaction and is none zero even 
for the parallel configuration of Fl and F2. 

Further over the last decade , the physical charac- 
teristics of the purity and entanglement of quantum- 
mechanical states has been recognized as a central re- 
sources in various aspect of quantum computation [14] 
and spintronics. Since for spintronics operation spin- 
orbit interaction is an integral part which entangles the 
orbital and spin degrees of freedom. This entanglement 
affects the coherence of quantum state in a nontrivial 
way. Therefore it is important to know how the purity 
and entanglement are affected by the spin-orbit interac- 
tion. Motivated by this we study in this paper, How does 
von-Neumann entropy (which is a measure of entangle- 
ment and purity of state) depends on various parameters 
of the problem? As we will see that this is also related 
with the torque arising due to SO interaction. Therfore 
seemingly different phenomena are related at fundamen- 
tal level. Our theoretical formulation provides a unified 
way to addresses all these question in one setting. 

In light of above discussion we develop a theory to 
study this effect quantitatively and qualitatively. Since 
the state of ferromagnet is not a coherent superposition of 
up and down electors rather it is a mixed state. Hence a 
proper scattering theory should take this mixed character 
of ferromagnetic state into account [15] . This is done con- 
veniently in density matrix formulation. To this end let 
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us assume that an electrons in Fl and F2 occupies states 
\n, a) = \n) ® \a) where n is channel index due to trans- 
verse confinement along y direction and a determines the 
spin state of the channels and can take tow values f and 
J. respectively. Let an electron be incident from Fl onto 
2DEG from a state tpi n = \n,a) while the final sate in F2 
is given as tp f = T\n,a) = £) m)j8 |m,/?)(m,/3|T|n,a) = 
Em^mnl™!^}- Here we have used the completeness 
relation J2 m l TO > P)( m i P\ = 1 ar± d ^mn i s transmission 
coefficient from state \n,a) to state m, 0). . The final 
state density matrix is given by 

Pf a = ^l)(^f\= E Tg a n Tlr\mp)(h\, (2) 

where N is normalization factor ensuring that Tr(p) = 1. 
Since we are interested in the spin degrees of freedom 
only hence we should take trace over orbital degree of 
freedom in eq. (2), which leads to the following equation 
for spin density matrix, i.e., 

m 0-y 

Since the incident electrons are also originating in 
mixed state so the corresponding density matrix for the 
incident electron is p,„ = n a \a) + — a), where n a 
and ri- a is the number of up electrons and down elec- 
trons respectively. The polarization P of ferromagnet is 
related to these as P = (n a — ri- a )/(n a + n_ a ). The 
fact that incident state is a mixed is incorporated into fi- 
nal state density matrix leading to the final state density 
matrix as, 

Pf=n a p] a + n_ a p n f - a (4) 

Once the final sate density matrix is known one can 
calculate the different component of transported polar- 
ization by taking respective trace of Pauli matrices i.e., 
ptr __ Tr(pfai), where i—x, y or z. 

Von-Neumann entropy for the final state density ma- 
trix is defined as, 

S v = -Tr(p f ln(p f )) = -J2 AjZn(Aj), (5) 

i 

where A^ are the eigenvalues of the density matrix pf. 
From the definition it is clear that extremal values of S v 
are and 1 for pure and completely mixed state respec- 
tively. 

To obtain quantitative results we perform numerical 
simulation for a two dimensional lattice of 100 x 100 sites. 
We model the conductor on a square tight binding lattice 
with lattice spacing a and we use the corresponding tight 
binding model including spin orbit interaction given by 
eq.(l) [17]. For the calculation of spin resolved trans- 
mission coefficient , we use the recursive green function 



method. Details of this can be found in Ref. [17], [18]. 
We fix the Fermi energy Ep — l.Ot above the bottom of 
band where t is tight binding hopping parameter we take 
t—1 as the unit of energy. The tight binding dimension- 
less SO parameter a = X/t * a where a is lattice spacing. 
The numerical result presented takes the quantum effect 
and multiple scattering into account. For the model of 
disorder we take Anderson model, where on-site energies 
are distributed randomly within [-U/2, U/2], where U is 
the width of distribution. 

Fig. 1 show the component of the transported polariza- 
tion i.e., Pf , Py T and Pf respectively for ballistic sys- 
tem as a function of absolute angle 9 of the magnetization 
direction of Fl and F2. We remind that since the system 
is two dimensional a natural frame is defined by plane 
of 2DEG and a normal to the plane in right hand sens. 
In this frame we rotate the magnetization of Fl and F2 
along the direction (0, — sin(0), cos#) while keeping them 
always parallel i.e. Pi = P2=Po(0, — sin(0), cos(0)). 
This implies that to begin with there is no component 
of polarization along x i.e. P x = and it remains so 
in absence of spin-orbit interaction. This is confirmed 
by numerical calculation where in Fig. 1 we have shown 
the results for zero SO interaction and as expected P x 
remains zero while P y and P z shows variation as — sin(#) 
and cos(#) as expected. Further we see from Fig. 1 
that Pf ^ for non zero strength of spin orbit inter- 
action as well the other two components along y and z 
axis. This confirms our prediction made in introduction 
that due to SO interaction, transported spin polarization 
will develop a component perpendicular to the magneti- 
zation direction of F2. Hence the torque acting on F2 
will be non zero even though Fl and F2 are parallel. 
This is seen clearly in Fig. 2 where torque is plotted 
as a function of absolute angle 9 for ballistic(C/ = 0.0) 
and diffusive([/ = 0.5t and U = l.Ot) case respectively. 
Dimensionless torque is defines as Torque — ^ j/| e |) I 
P 2 x P tr |= P VTPf 2 + (P* r cos(0) + Pf sin(0) 2 ). We 
notice that for ballistic case the torque is zero when Pi 
and P2 are parallel to y axis. This is so because injected 
electrons in this case are approximately an eigen state of 
RSO interaction given by eq.(l). While for the diffusive 
case torque is non zero for all angles 9, since in diffusive 
case due to the scattering direction of Br(1c) is random- 
ized hence probability of finding electron in eigenstate of 
SO interaction is reduced. 

Left Panel of Fig. 3 shows magnitude of net trans- 
ported polarization i.e., P tr = ^[Pf 2 + P* r 2 + Pf 2 ) 
for ballistic system with varying SO interaction strength 
a . We notice that for = 90, P tr is almost equal to the 
value of polarization in the absence of spin-orbit interac- 
tion , shown as straight line. This is consistent with the 
results of Fig. 2 where torque vanishes for 9 = 90. Fig. 
3 (right panel) shows entropy for ballistic system as a 
function of absolute angle 9. First we notice that S v < 1 
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for a — 0.0 implying that the states of Ferromagnetic 
leads Fl and F2 are a mixed state. Therefore the use of 
density matrix scattering theory is essential and the only 
way to treat these mixed state situation. Further as is 
seen entropy is minimum when P r is maximum reflecting 
the more pure character of state and is maximum when 
P tr is minimum reflecting the more mixed character of 
the state. Moreover the entropy can be less compared to 
the vale for zero SO coupling for certain values of 9 while 
corresponding P tr is large compared to the value for zero 
spin orbit. Implying that SO interaction has a polarizing 
effect and can in principal create a pure state from mixed 
state [19]. This is important since in general any scat- 
tering process is reverse that is pure states evolve into 
mixed sates and one looses the associated quantum char- 
acter of states which is rather important for spintronics 
and quantum computation. 

In Fig. 4 we study Torque and Entropy as a func- 
tion of exchange splitting. We see that torque remains 
zero when polarization of ferromagnet points along y axis 
(9 = 7r/2). This is in agreement with the results shown 
in Fig. 2 and arises since the injected electrons are ap- 
proximately in the eigen state of SO interaction. While 
the torque along z axis (9 = 0.0) first increases with the 
increase of A and then drops to zero at A/t = 2.0. At 
the same same point entropy also goes to zero implying 
that the injected state is pure. Fig. 5 confirms that at 
A/t = 2.0 state becomes pure since polarization goes to 
one. In other-words it is not entangled with the orbital 
degrees of freedom. Since the state is decoupled from 
orbital degrees of freedom hence the torque arising due 
to spin-orbit interaction which couples spin and orbital 
degrees of freedom goes to zero. Therefore for pure state 
the torque arising due to SO interaction will be zero irre- 
spective of along which direction it is pointing. Further 
this also implies that torque is anisotropic and depends 
on absolute angle. 

Fig. 6 and Fig. 7 shows Torque ,P tr , and Entropy 
respectively along Z axis as function of SO interaction a. 
We see that all three quantities shows oscillation as a is 
varied and the period of oscillation 8a — 0.03, which cor- 
responds to precession angel of tt over length 100 a. This 
length is the length of our systems. Hence these oscilla- 
tion are due to spin precession caused by SO interaction 
over the system length. 

In conclusion we have shown that in collinear spin valve 
structure a non zero torque arises due to the presence 
of SO interaction. We have related this with the Von- 
Neumann entropy and studied their behaviors as function 
of different parameters. We have formulated a spin den- 
sity matrix scattering theory which allows us to study all 
these seemingly different phenomena under one setting. 




FIG. 1. Different comopnent of transported spin polariza- 
tion as a function of absolute angle 6 for different values of 
dimensionless SO interaction strength a. Other parameters 
are i?F = 1.0,and exchange splitting(A), in Ferromagnet Fl 
and F2 is given as 2 * A/ Ep — 1.5. 
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FIG. 2. Dimensionless torque (defined in text) as a func- 
tion of absolute angle 6 for different values of dimensionless 
spin-orbit coupling a. Other parameters are same as for Fig. 
2. 




FIG. 3. Transported polarization (left panel) and 
von-Neumann entropy as a function of angle 9 for different 
values of dimensionless spin-orbit coupling a. Other parame- 
ters are same as for Fig. 2. 
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FIG. 4. Dimensionless torque (defined in text) and Entropy 
along y and z axis as function of exchange splitting A. Other 
parameters are same as for Fig. 1. 



FIG. 7. Transported polarization (left panel) and 
von-Neumann entropy as a function of of dimensionless 
spin-orbit coupling a. Other parameters are same as for Fig. 
2. 
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FIG. 5. Polarization along y and z axis as function of ex- 
change splitting A. Other parameters are same as for Fig. 
1. 
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FIG. 6. Dimensionless torque (defined in text) along z axis 
as a function of dimensionless spin-orbit coupling a. Other 
parameters are same as for Fig. 1. 
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